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summary 

A  method  of  solving  minimal  cost  network  flow  problems  Is 
described.  Tills  method  begins  with  any  circulation,  feasible 
or  not,  ana  an  arbitrary  pricing  vector.  A  labeling  procedure 
Is  then  usea  10  adjust  •out-of-kilter*  arcs,  that  is,  arcs 
tnat  fail  to  satisfy  the  optimality  properties.  It  is  shown 
that  the  mot  ou  terminates  in  a  finite  number  of  steps,  ana 
that  in  so  doing,  the  status  of  no  arc  of  the  network  (as 
measured  by  certain  •filter  numbers*)  is  worsened  at  any  step. 
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AM  0 UT—O F— X I LTE R  METH OD  FOR  MINIMAL 
COST  FLO,;  PROBLEMS1 
D.  R.  Fulkerson 

1.  Introduction  In  thin  paper  a  method  of  solving 
minimal  cost  network  flow  problems  is  described  and  shown  to 
have  a  certain  desirable  monotone  behavior.  The  method  begins 
with  an  arbitrary  flow,  feasible  or  not,  together  with  an 
arbitrary  pricing  vector,  and  then  uses  a  labeling  procedure 
to  adjust  an  arc  of  the  network  that  falls  to  satisfy  the 
appropriate  optimality  properties. 

To  oresent  the  oasic  notions  underlying  the  method,  let 
ur  consider,  for  a  moment,  a  general  linear  program  of  the 
form 


(1.1) 

1 i 

(1.2) 

M  <  *, 

(1.3) 

minimis 

■fere  the 

X  -  (Xp 

au’  V 
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x  Is  few. 

1  Tills  res 
System. 
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X,  -  b 
d  X 
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vector  v  =  (ir,  ,  v  )  such  that  the  implications 

1  m 

(1.4)  0j  +  >  0 

m 

(1.5)  c,  +  £  Tr.a  .  <  0 

u  i^l  1 

hold  for  all  j.  Then  It  follows  that  x  is  a  minimizing  solution, 
and  thus  (1.4),  (1.5)  might  be  termed  optimality  properties. 

For  a  given  x  satisfying  (l.l)  and  for  any  i r,  the  following 
case  classification  for  the  j— th  component  of  the  program  is 
exclusive  and  exhaustive: 

°j  +  '  Vlj  >  °>  *3  '  l3 

°j  +  V  Vij  “  l3  *  *3  *  us 
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J 
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73  i^i  J 

<  0, 

XJ 

<  u, 
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*  i 
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*  0, 

x. 
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If  all  comno.nents  are  in  one  of  the  states  a,  £,  7,  then  x  is 
feasible  and  optical.  /e  call  these  the  "in— kilter’*  states,  the 
others  "out— of-kilter"  states.  The  algorithm  to  be  presented  for 
network  flow  problems  concentrates  on  a  particular  out-of-k liter 
co^oonent ,  ana  gradually  put3  it  In  kilter.  It  does  this  in  such 
a  v/ay  that  all  in— kilter  components  stay  in  kilter,  whereas  any 
other  out-of-kilter  comoonent  either  improves  or  stays  the 
same,  in  a  sense  made  precise  in  section  2. 

Section  2  provides  a  description  of  the  special  class  of 
linear  programs  to  which  the  method  applies,  together  with  some 
preliminary  discussion.  assume  that  the  given  data  for  the 
program  are  integers  (or,  equivalently,  rational).  Then  the 
algorithm,  presented  in  section  2,  works  with  integers  throughout 
A  oroof  that  the  algorithm  terminates  in  a  finite  number  of  steps 
and  ti  at  in  so  doing  it  possesses  the  monotone  property  roughly 
described  above,  In  sketched  in  section  4.  Some  comparisons 
with  other  methods  for  solving  minimal  cost  flow  nroblems  are 
made  in  section  5» 

For  the  particular  class  o f  programs  being  considered,  the 
assumption  that  the  initial  x  satisfies  (i.l)  is  unimportant, 
since  such  an  x  is  Immediately  available,  e.g.,  x  =  0.  But 
starting  with  a  good  guess  for  x  and  v  will  decrease  comoutution 
time.  One  situation  for  which  the  nresent  algorithm  is  par¬ 
ticularly  appropriate  would  be  in  solving  a  sequence  of  flow 
problems,  where  each  problem  of  the  sequence  differs  only 
slightly  from  Its  predecessor.  Then  the  old  optimal  x  and  i r 
could  be  used  to  initiate  the  comoutation  for  the  new  Qroblem. 
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We  should  like  to  express  our  appreciation  to  G.  B. 
Dantzig,  whose  criticism  of  an  earlier  version  of  this  paper  in 
which  the  intltial  x  was  assumed  feasible,  led  U3  to  reconsider 
the  problem  from  the  standpoint  of  infeasible  x. 

2.  Notation,  definitions,  and  problem  description. 

We  suppose  given  a  network  consisting  of  nodes  1,  2,  . ..,  n 
together  with  directed  arcs  ij  (from  node  i  to  node  J ) •  Each 
arc  1J  has  associated  with  it  three  integers:  (the  arc 

lower  bound) ,  u^  (the  arc  upper  bound  or  capacity) ,  and 
c^  <  (the  arc  cost) ,  with  0  £  ,  £  u1 

It  is  convenient  to  describe  the  problem  in  terms  of  cir¬ 
culations  fl2] ,  rather  than  flows  from  sources  to  sinks 
i7, 8,  a.  By  a  circulation  we  shall  mean  a  nonnegative 
integral  vector  x  *  (x.  ,),  one  component  for  each  arc  1J,  that 

J 

satisfies  the  conservation  equations 

( —  • 1 )  —  xj^)  =  0  (i  =  1#  •••»  n) . 

If  the  circulation  x  also  satisfies 

(2.2)  (all  arcs  ij), 

we  call  x  a  feasible  circulation.  We  shall  refer  to  a  parti¬ 
cular  component  of  a  circulation  as  the  arc  flow  or 
the  flow  in  arc  1J. 

A  feasible  circulation  x  that  minimizes  the  cost  form 


(2.3) 


7. 

1.1 


CiJXU 


over  all  feasible  circulations  is  optimal.  The  problem  we 
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are  considering  is  that  of  constructing  an  optimal  circulation. 
Of  course  feasible  circulations  may  not  exist,  in  which  case 
te  .-/ant  to  discover  this  fact.  It  is  known  [12]]  that  a  neces¬ 
sary  and  sufficient  condition  for  the  existence  of  a  feasible 
circulation  is  that  the  inequalities 


<2. 4) 


y  U11  2  1 

1J  1  in  L  J1 

3  In  I 


i  in  L 
J  in  Z 


hold  for  all  subsets  L  of  nodes.  Here  T  denotes  the  complement 
of  L.  The  conditions  (2.4)  are  easily  shown  to  be  nece33ary; 
their  sufficiency  can  be  proved  in  various  ways,  for  example, 
using  the  max  flow  min  cut  theorem  [j,  8]  or  the  supply- 
demand  theorem  [ll[] . 


Let  7 t  -  (7^)  be  a  vector  of  integers,  one  component  for 
each  node  i.  «’e  call  v  a  pricing  vector,  and  refer  to  its 
components  as  node  prices .  Optimality  properties  for  the 
problem  are  that  the  implications 


(2.5)  +  TTi  -  TT1  >  0  - »X^  ^  ij, 


(2.6) 


ij 


m  m  < 0  — *  =  ui 


’no Id  for  all  arcs  ij.  That  is,  if  x*is  a  feasible  circulation, 
and  if  there  is  a  pricing  vector  tt  such  that  (2.^),  (2.6) 
hold,  then  x  is  optimal.  We  shall  shorten  the  notation  by 
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(2.V)  C1?  C1;  -4  TTj  -  T,. 

F'r  a  given  circulation  x  and  pi’icin?  ector  t,  an  arc. 
1,1  is  in  .hist  one  of  ♦’he  following  states: 
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the 

arc  'liter  numbers  corresponding  to  out— of— !:1  Iter  states  are 
listed  ielow: 


(a.  ) 


/  —  y 

Xi,1 


4  i 
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(M 

j. 

*i:  -  *1.1 

(7  ) 

'i.-W;  -  ui.  > 

(a  ) 

7 1  d x i 

(3  .) 

‘“At  ~~  ’  *  A  t 

c. 

1. 

(>?) 

V  —  It 

j  ui  « 

Thu;,  out-of— '  liter  ar  :s  • .  a  ve  cosi  tive  Miter  numbers .  The 
'il’.rr  numbers  for  states  a,,  ,  {3.,,  *y0  measure  ln^easlblll  :y 

fo.  the  ;:ro  flow  ,  w:.lle  the  •'liter  numbers  for  states 

v-‘ 

Y-,  .  a0  are  a  measure  of  the  decree  to  which  tne  optimality 


r* »p4  *  0 

Ac  C,\ 

(?.6) 

"all 

to  be  satisfied. 

Th-  a 

Igcrj  th: 

■  .3  r  .ut 

3  J  1  n 

the  following  section  lias  the 

pr^oertj-  t 

hut  all 

arc  . 

1  Iter 

a. mberu  are  monotone  nonincreasing 

t  ro ■.*;?  f,jr 

*  y  a  •> 

,1  l  V  /  » 

nputo c 

1  0  f .  • 

however,  steps  "can  occur  that 

otht'.re  no  liter  nur  her .  a  ia  i.r.is  complicates  the  proof  o1' 
termination  somew  ;fct. 

need  a  few  other  notions  before  stating  the  algorithm, 
the  main  one  being  that  of  a  path  from  some  node  to  another  In 
a  network.  Let  1^,  i2,  lin  be  a  sequence  of  distinct 

nodes  of  a  network  such  that  either  1,1.,.  or  1.  ,  1. 
is  an  arc,  k  =  1,  m  —  1.  Picking  out,  for  each  k,  one  of 

these  two  possibilities,  we  call  the  resulting  sequence  of 
noies  and  arcs  a  path  from  1^  to  iJn.  Arcs  that  belong 

to  the  path  are  forward  arcs  of  the  path;  arcs  that 


belong  to  the  path  are  reverse  arcs  of  the  path.  If  we  alter 
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the  definition  of  a  oath  b;  stimilatirg  tna*  3-  1  ,  we 

i  m 

call  tne  resulting  sequence  of  nodes  and  arc:  a  cycle . 


An  out-s>f "filter  algorithm.  The  algorithm  of  this 
section  uses  a  modified  labeling  procedure  [B,  a]  as  its 
basic  routine.  In  general,  the  labeling  procedure  is  a 

search  for  u  pat!:  (having  certain  desired  properties)  fro::, 
some  node  to  another.  Ve  3 tart  labeling  fro*:,  a  given  node, 
cal ’ea  the  origin .  attempting  to  reach  some  other  given  node, 
called  the  terminal .  To  1  iltiate  the  modified  procedure,  we 
assign  tne  laoel  [p ,  to  the  origin:  the  following  labeling 

rules  :r«  then  applied: 

(  :.l)  If  node  1  Is  labeled  Q-.-,  e.j,  node-  *  ii 

uni '-.be  led ,  and  if  1.'  Is  ar  crn  suci  t  .at  either 
(a)  h  |  y  »  * «  *  \  t  j  ( » 

:  X.  .  X.  , 

V  •  •  t  '  1  '  Si  *  '  ‘  ^  '  N  u  £  •  • 

then  node  r-.ceives  tne  atel  J1  ‘ ,  j  ,  where 

-  ::.in(r,,  -  x<  .)  I  ca.>e  (a),  ,  -  mir.^,  u±  j  -  X^) 

in  case  (  ). 


1“  ode  I  1  l’teled  t  «j  ,  node  '  Is  unlabeled, 

an]  I*'  '1  is  ar  or?  eh  the  t  either 


(a)  ' 


*i  - 


V.' 


T 
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(h)  cn  <  0,  x1jL  >  uJ4 


then  node*  j  receives  the  label  [l  ej,  where 
e1  -  *ain(  -  1^)  in  ca3e  (?  ) , 

e1  -  min(c^,  x  —  u^)  In  case  (b). 


Here  Is  a  circulation  and  it  a  pricing  vector. 

Ti  e  labeling  procedure  terminates  in  one  of  two  ways, 
called  breakthrough  ar»d  nonbreak through ,  respectively:  either 
the  terminal  receives  a  label,  or  no  more  labels  can  be 
assigned  and  the  terminal  nas  not  been  labeled. 

If  breakthrough  occurs,  a  path  from  origin  to  terminal 
can  be  located  by  bacictracxing  from  the  terminal,  using  the 
first  ■‘.ember:  of  the  label  pairs.  If,  In  this  backtracking, 

&  node  J  is  reached  that  carries  the  label  JjL+,  c  J,  then  1 J 
is  a  forward  arc  of  the  path  from  origin  to  terminal;  If  .1  is 
labeled  [i“",  €  ^  >  then  J1  13  a  reverse  arc  of  the  path.  Thus 

t 

forward  arcs  of  the  path  satisfy  either  (3* la)  or  (3»lb)» 
whereas  reverse  arcs  of  the  path  satisfy  (5.2a)  or  (3.2b). 

If  nonbreakthrough  results,  we  let  L  and  T  denote  the 

i 

sets  of  labeled  and  unlabeled  nodes  respectively ,  and  define 
two  subsets  of  arcs : 


(3.3) 


xu  £  ui<] 
hi  ^  'n}- 


(■5.*) 
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We  then  define 

(3.5)  6,  -  nin  (c.  ) 

1  u  in  a2 

(3.6)  50  =  min  (-  c  M ) 

J1  In  J1 

(3.7)  5  -  mlnftj,  r>2). 

Here  (l  =  1,  ?.)  Is  a  positive  Integer  or  00  according  ao 

Is  nonempty  or  empty.  * 

The  complete  algorithm  now  runs  as  follows.  Start  the 

computation  with  any  circulation  x  and  any  pricing  vector  v. 

Next  locate  an  out-of— !: liter  arc  st  and  go  on  to  the  appro¬ 
priate  case  below: 

1 

(a,)  c3t  >  0,  xst  <  fst.  The  origin  for  labeling 
Is  t,  the  terminal  s.  If  breakthrough  re¬ 
sults,  add  c  -  mln(e„ ,  -  xct)  to  the 

flow  In  all  forward  arc 3  of  the  path  from 
t  to  3,  subtract  t  from  the  flow  in  all 
reverse  arcs,  and  add  c  to  xst«  If  non- 
breakthrough  results,  add  5  defined  in  (3.7) 
to  all  iri  for  1  In  L. 

(9j)  or  (Tj)  est  =  0.  *3t  <  <st.  or  ?5t  <  0. 
xs;  <  ust-  Same  as  (Oj),  except 
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t  -  ">m(V  ust  “  xst)- 

(a2)  or  (e2)  cgt  >  0,  xst  >  igt,  or  cgt  =  0, 

x  .  /  u  .  The  origin  for  labeling  in  s, 

T*  u  l 

the  terminal  t.  If  breakthrough  results, 

acd  .  min(c.  ,  x  .  —  I  )  to  the  flow  In 

all  forward  arcs  of  the  path  from  a  to  t , 

subtract  c  from  the  flow  In  all  reverse 

arc:'  and  subtract  t  from  x  If  nonbreak— 

s  t 

through  result  ,  add  b  defined  In  (5.7)  to 
all  v ,  for  i  in  Th 

(yr)  cgr  \  0,  xGt  /  ur;t.  Same  as  (a^)  or  (0?), 
except  *:  -  n;in(t,_,  x  .  —  u  .  ) . 

The  labeling  process  is  reseated  for  the  arc  at  until 
eit  er  st  is  in  kilter,  or  until  a  nonbreakthrough  occurs 

I 

for  w:  lc  the  node  urlce  change  o  -  oo  .  In  the  latter  case, 
stop.  (There  is  no  feasible  circulation).  In  the  former 
case,  locate  another  out— of— .cl Iter  arc  and  continue. 

4.  Termination  and  the  monotone  property.  Suppose  that 
arc  st  is  out  of  kilter,  say  in  state  a^.  The  origin  for 
labeling  lc  t,  tr.e  terminal  s.  The  arc  st  cannot  be  used  to 
label  s  directly,  since  neither  (5.  fa)  nor  (.3.?b)  is  applicable. 
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Consequently,  if  breakthrough  occurs,  the  resulting  path  from 
t  to  s,  together  with  the  arc  st,  is  a  cycle.  Then  tie  flow 
changes  that  are  made  on  arcs  of  this  cycle  again  yield  a 
circulation.  Moreover,  the  labeling  rules  have  been  selected 
in  such  a  way  that  Icilter  numbers  for  arcs  of  this  cycle  do 
not  increase,  and  at  least  one,  namely,  for  arc  st.  decreases. 
Kilter  numbers  for  arcs  not  in  the  cycle  o'1  course  don't 
change . 

Similar  remarks  apply  if  st  is  In  one  of  the  other  out— 
of— Miter  states. 

he  :ul’ imari  :e  toe  possible  effects  of  a  breakthrough  on  an 
arc  i;  In  t  :c  diagrui.  :>elo  t ,  which,  snows  the  state  transitions 
that  :.;ay  occur  following  breakthrough.  If  a  transition  is 
possible,  the  number  recorded  b  tide  the  corresponding  urrew 
represents  e  change  in  kilter  maker.  (The  subscripts  i,: 
arc  omitted  in  the  uiagran.) 


ft  ,->r*  —  -n  ft  O  n  — : 


Brea'  through  Diagram 
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Verification  of  the  breakthrough  diagram  is  straightforward 

For  example,  suppose  arc  ij  is  in  state  a^,  with  /  0, 

xi  i  ^  ^ij»  anri  filter  number  c^^x^  —  X^j)  0*  If  IJ  is  not 

an  arc  of  the  cycle  of  flow  changes,  then  ij  remains  in  state 

with  zero  change  in  kilter  number.  If  the  flow  In  arc  IJ 

has  changed  as  a  result  of  the  breakthrough,  then  either  ij  is 

the  arc  at  or,  by  the  labeling  rules  (3*l)»  (3*2),  IJ  Is  a 

reverse  arc  of  the  oath  from  origin  to  terminal.  Specifically, 

i  was  labeled  from  J  using  (3.2a).  In  either  case,  x,  .  de— 

ij 

. 

creases  by  the  positive  integer  £  £  xi  j  ”  thc  ne*  8tate 

for  ij  is  Op  or  a,  and  hence  the  kilter  number  for  ij  has 
decreased  by  .  e  )  0.  The  rest  of  the  diagram  may  be 
verified  similarly. 

The  state  transitions  and  changes  in  kilter  number  that 
may  occur  following  a  nonbreakthrough  with  6  a>  are  Indicated 
in  he  following  diagram.  (Again  the  subscripts  Ij  are 
omitted. ) 
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Again  we  omit  a  detailed  verification,  but  consider,  for 
exanole,  an  arc  i*  in  state  7^,  so  that  ^  <  0,  x^j  <  u^, 
having  kilter  number  —  ulj)  ^  ®  before  the  node  price 

change  is  made.  If  both  i  and  J  are  in  L  or  both  in  t,  then 


i  remains  the  same  after  the  node  crlce  change,  and  conse¬ 
quently  ij  stays  in  state  7^  with  no  change  in  kilter  number. 

Are  cannot  have  i  in  L,  j  In  t  (labeling  rule  (3.1b)),  and  hence 
the  remaining  possibility  is  i  in  L,  j  In  L.  Then  A  is 
increased  by  6  )  0.  Consequently  the  arc  1,1  either  remains  in 
state  Yj*  (if  6  <  -  Cj.),  goes  into  state  3  (if  fi  =  and 

2  ijj)*  into  state  P1  (if  6  =  -F^.  and  <  ljj)»  or  into 
state  a1  (if  5  )  and  xi «  ^  i^,)#  and  the  corresponding 


El 
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chances  in  I'llter  number  are  respectively  c;(x,  —  u,  .)  <(  0, 

3(x.  ,  -  uj  ,)  <  0,  £11  -  xt  <  -  6(x1  ,  -  u.  .)  <0, 

'  ~  Ki  '  “  ^1  ,(x^  ~  u ,  J  ^  0.  (The  remaining  logical  po33ibi- 

11  ty  s  /  -e,,,  ,  2  •  cannot  occur,  since  if  x^,  /_t±y  then 

ia  is  In  61-  defined  bv  (y.4),  and  lienee  f  K  -c,,.) 

<-  ~  t  J 

It  follows  fro  i  trie  breakthrough  and  nonbreakthrough 

ilarramr  that  kilter  numbers  are  monotone  nonincreasing 

througnout  the  comoutat ion .  Moreover,  if  breakthrough  occurs, 

at  least  one  arc  ’'liter  number  decreases.  Thus  to  prove  that 

the  algorithm  terminates,  it  suffices  to  show  that  an  infinite 

sequence  of  successive  nonbreakthroughs,  each  with  5  <(  oo  ,  is 

impossible.  To  show  this,  let  us  suppose  that  a  labeling 

resulting  in  nonbreakthrough  with  r  f  oo  has  occurred,  and 

let  L,  Z  denote  the  labeled  and  unlabeled  3ets  of  nodes.  After 

changing  node  prices,  the  new  cT  vector,  which  we  denote  by 
_  ? 

c  ,  has  components  given  In  terms  of  the  old  by 


r 


c  —  * 

A  ' 

X. 

f’or  1  in  L, 

J  in 

z 

(*.1) 

"  11 

c .  fc 

'  \  1 

1  1 , 

^or  i  In  Z, 

f  in 

L 

\ 

1 

*1,1 

otherwise . 

If  the  arc  st  is  .-till  oui  of  liter,  then  the  origin  is  the 
same  f'or  the  next  labeling,  and  jt  follows  from  (4.1)  and  the 
luoelJng  rules  that  every  node  o p  L  may  again  be  labeled. 
TIP'S  If  tic  new  labeling  again  results  in  nonbreakthrough 
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'  '  *  r 

with  labeled  set  L  ,  we  have  L  c  L  .  Let  Oj ,  u0  denote  the 

t  w  I 

new  jets  defined  in  terms  of  L  ,  c  (and  x)  by  (3.5),  (^.4), 
and  suppose  L  =  L  .  T.oen,  from  (4.l)  we  have  Ci1  ^  , 

Cl ^  c  Q^,  and  at  least  one  of  these  Inclusions  is  proper  l;, 
(3»5) >  (3*6),  (3*7) •  Hence  the  new  labeling  either  assigns  a 

label  to  at  least  one  more  node,  or  failing  this,  an  aro  Is 
removed  from  one  of  the  sets  L  or  It  follows  that,  after 


finitely  many  nonbreakthroughs  with  6  <  00 ,  we  either  get  the 

arc  st  in  hi] ter,  obtain  a  breakthrough,  or  obtain  a  nonbreah— 
through  with  5  =  00  . 

If  a  nonbreakthrough  with  6  00  occurs,  then  there  is 

no  feasible  circulation.  For  if  *  -  ao  ,  it  follows  from  (3*?)> 

(3.4)  and  the  labeling  rules  (3-l),  (3*2)  that  x^  ui  ^  ^or> 

i  In  L,  1  In  Z,  and  x..  <(  i  . ,  for  1  in  L,  j  in  XL  Moreover1, 

1-1  ,M 


for  the  arc  st,  either  t  is  in  L,  1  in  L  'with  x  ,  <  £,  .  ,  or  s  is 

o  1/  o  U 

in  L,  t  in  Z  with  x  .  )  u  .  (Thii  Is  immediate  for  cases 

o  L  u  l 

a.  ,  £, ,  7,,  of  the  alrorltnm,  and  follows  from  (y.3)  and 

the  assumption  h  ■-  00  for  case  o0,  from  (p.4)  and  the  assumption 
r  -  x)  for  case  7^ .  )  Hence,  summing  tre  equations  (2.1)  over 
1  in  L  and  noting  cancellations,  we  obtain  in  all  cases 


0 


V 

/ 

i  in  L 
J  in  I 


xn> 


i  In  L 
J  In  L 


But  this  violates  tne  feasioil^t'  condition  (2.4).  Thus 
6  ‘To  implies  there  is  no  feasible  circulation. 
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To  sun  up,  the  algorithm  terminates  after  finitely  many 
applications  of  the  labeling  procedure,  either  with  all  arcs 
in  kilter  (In  which  case  the  feasible  circulation  is  optimal), 
or  with  the  conclusion  that  there  is  no  feasible  circulation. 
Moreover,  all  arc  kilter  numbers  are  monotone  nonincreasing 
throughout  the  computation. 

It  is  worthwhile  to  note  the  simplification  that  occurs 
if  the  method  of  the  preceding  section  i8  initiated  with  a 
feasible  circulation.  The  states  a^,  y ^  are  then  empty 

to  begin  with,  and  consequently  remain  empty  throughout  the 
computation.  Hence  at  each  nonbreakthrough  (as  well  as  each 
breakthrough),  the  kilter  number  for  at  least  one  arc,  namely 
st,  decreases  by  a  positive  integer.  In  many  minimal  cost  flow 
problems,  a  starting  feasible  circulation  is  readily  at  hand. 
For  example,  in  the  Hitchcock  problem  [l,  2,  3,  4]  or  the 
assignment  problem  Qj,  6,  10] ,  such  is  the  case. 

5»  Comparison  with  other  methods.  The  method  of  section 

3  Is  a  generalization  of  the  method  of  [9]  for  solving  minimal 

cost  flow  problems,  which  in  Itself  generalizes  the  methods 

of  [5»  6,  8,  10]  for  solving  Hitchcock  and  assignment  problems. 

In  [Q]  the  fundamental  problem  was  that  of  finding  a  maximal 

feasible  flow  from  source  node  1  to  3lnk  node  n  that  minimizes 

cost  over  all  such  flows.  (Also  the  lower  bounds  were  assumed 

zero  on  all  arcs.  Thin  Is  not  really  a  restriction,  since  a 

change  of  variables  will  accomplish  this,  if  desired.)  If  we 

add  to  the  network  the  special  arc  nl  with  l  ,  *  0.  u  , 

nl  '  nl 


=  U, 
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c  ^  =  —  C  (U  and  C  lcrge),  and  consider  feasible  circulation' 
In  the  enlarged  network,  then  the  method  of  section  3  is 
aoplicable  to  such  problems.  Or  if  it  is  desired  to  find  an 

optimal  flow  from  1  to  n  of  given  value  v  -  7(x,  .  —  x.,)  in 

J  J  J 


the  original  network,  we  can  add  the  arc  nl  with  i  ,  =  u  ,  -  v, 

ni  nl 

cnl  =  order  t0  the  problem  in  circulation  form. 

The  method  of  [9]  begins  with  the  zero  flow  from  source  1 
to  sink  n  (which  satisfies  the  bounds  on  arc  flows  because 
lOv  r  bounds  are  zero),  and  all  node  prices  zero.  It  was  also 
assumed  that  the  given  arc  costs  are  nonnegative.  Equivalently, 
if  we  take  £nl  *  0,  un^  =  U,  c  ^  =  —  C  and  begin  the  algorithm 
of  section  3  with  the  zero  circulation  and  all  node  Drices 
zero,  then  the  special  arc  nl  is  the  only  out-of— kilter  arc 
(it  is  in  state  7^),  and  hence  it  remains  tie  only  out— of-k liter 
arc  throughout  the  computation.  Then  the  method  of  section  ;> 
reduces  to  that  of  [9] . 

It  is  also  informative  to  note  some  of  the  major  contrasts 
between  this  method  and  the  simple/  method  |_4j  for  solving  such 
problems.  First  of  all,  the  siiaolex  method  would  be  done  in 
two  phases,  the  first  phase  being  a  search  for  a  feasible 
circulation,  the  second  for  an  optimal  circulation.  (Through¬ 
out  both  of  these  phases,  the  simplex  method  would  work  with 
basic  solutions,  a  concept  that  play3  no  role  in  this  method.) 
Here  we  hsve  combined  the  two  phases.  Ignoring  this  difference, 


however,  and  assuming  that  both  methods  start  with  a  feasible 
circulation,  the  main  contrast,  apart  from  mechanics  of  operation. 
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appears  to  lie  in  the  fact  that,  for  the  simplex  method,  the 
kilter  numbers  are  not  nonotone.  For  example,  arcs  that  v/ere 
in  kilter  at  some  stage  of  the  simple::  computation  can  go  out 
of  kilter  at  later  stages. 
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